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SUMMARY 


The film-parameter equations are provided for four fluid-film lubrication regimes 
found in elliptical contacts. These regimes are isoviscous-rigid; viscous-rigid; 
elastohydrodynamic of low-elastic-modulus materials, or isoviscous-elastic; and 
elastohydrodynamic, or viscous-elastic. The influence or lack of influence of elastic 
and viscous effects is the factor that distinguishes these regimes. The film-parameter 
equations for the respective regimes come from earlier theoretical studies by the 
authors on elastohydrodynamic and hydrodynamic lubrication of elliptical conjunctions. 
These equations are restated and the results are presented as a map of the lubrication 
regimes, with film-thickness contours on a log-log grid of the viscosity and elasticity 
parameters for five values of the ellipticity parameter. These results present, for the 
first time, a complete theoretical film-parameter solution for elliptical contacts in the 
four lubrication regimes. 


INTRODUCTION 


Two solids, as shown in figure 1, having different radii of curvature in a pair of 
principal planes (x and y) passing through the contact between the solids, make con- 
tact at a single point under the condition of no applied load. When the two solids have 
a normal load applied to them, the point expands to an ellipse, with a being the semi- 
major axis and b the semiminor axis. It has been common to refer to elliptical con- 
tacts as point contacts. But, since this paper deals with loaded contacts, the term 
elliptical contact is used. For the special case where r^y = r^ and r-gy = rg x , the 
resulting contact is a circle rather than an ellipse. Furthermore, as either or 

rgy becomes much larger than the remaining radii of curvature, the elliptical contact 
approaches a rectangular contact. A rectangular contact exists when side leakage of 
the fluid can be ignored in comparison with the fluid flow in the rolling direction. 
Therefore, this paper is concerned with conjunctions between solids with contact re- 
gions ranging from circular to rectangular. 

The types of lubrication within the contact shown in figure 1 are influenced by two 
major physical effects: the elastic deformation of the solids under an applied load, and 
the increase in fluid viscosity with pressure. Therefore it is possible to have four re- 
gimes of fluid-film lubrication, depending on the magnitude of these effects and on 



their relative importance. These four regimes are defined as 

(1) Isoviscous-rigid: In this regime the elastic deformation of the surfaces is an 
insignificant part of the fluid film separating them, and the maximum pressure in the 
contact is too low to significantly increase fluid viscosity. This form of lubrication is 
typically used in circular-arc thrust bearing pads; industrial coating processes in 
which paint, emulsion, or protective coatings are applied to sheet or film materials 
passing between rollers; and very lightly loaded cylindrical roller bearings. 

(2) Viscous-rigid: In this regime the pressure within the contact is high enough so 
that there is a significant increase in the fluid viscosity within the contact. The defor- 
mation of the surfaces is an insignificant part of the fluid-film thickness. This form of 
lubrication is used on roller end-guide flanges, in contacts in moderately loaded cylin- 
drical and tapered rollers, and between some piston rings and cylinder liners . 

(3) Isoviscous-elastic: In this regime the elastic deformation of the solids is a 
significant part of the thickness of the fluid film separating them. The pressure within 
the contact is low and close to that obtained under dry frictionless conditions. This 
form of lubrication is used in seals, human joints, tires, and elastomeric-materigl 
machine elements. These applications are typified by solid materials with low elastic 
modulus. 

(4) Viscous-elastic: In this regime the elastic deformation of the solids is also a 
significant part of the thickness of the fluid film separating them, and the pressure 
within the contact is high enough to cause a significant increase in the fluid viscosity 
within the contact. This form of lubrication is used in ball and roller bearings, gears, 
and cams. 

Several authors (refs. 1 to 6) have contributed solutions for the film thickness in 
the four lubrication regimes, but their results have been confined largely to rectangular 
contacts. The essential difference between these contributions is the way the parame- 
ters were made dimensionless. In the present paper the film thickness is defined for 
the four fluid-film lubrication regimes just described for conjunctions ranging from 
circular to rectangular. The film -thickness equations for the respective regimes come 
from theoretical studies reported in references 7 to 9 on elastohydrodynamic and hydro- 
dynamic lubrication of elliptical conjunctions. The results presented herein are valid 
for isothermal fully flooded conjunctions. In addition to the film -thickness equations 
for the respective regimes, a map is presented of the lubrication regimes, with film- 
thickness contours on a log-log grid of the viscosity and elasticity parameters for five 
values of the ellipticity parameter. 
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SYMBOLS 


a 

b 

E 

E* 


F 

G 

g E 


g V 

H 

H 

h 

k 

Piv, as 
PO 

Of 

R 

r 

U 

u 

W 


x,y 

a 


semimajor axis of contact ellipse 
semiminor axis of contact ellipse 
modulus of elasticity 


2 



normal applied load 

dimensionless material parameter, E’/p^ v ag 

dimensionless elasticity parameter, W 8 ^ 3 /U 2 

o o 

dimensionless viscosity parameter, GW /U 


dimensionless film thickness, h/R 

A 


w \ 2 _ 


dimensionless film parameter, H( — ) = 

\V 


A 

2 2 d 3 
u ^0 R x 


film thickness 

ellipticity parameter, a/b 

asymptotic isovisco us pressure 


maximum Hertzian pressure 
reduced pressure parameter 


effective radius 
radius of curvature 

dimensionless speed parameter, utJq/E'R x 
mean surface velocity in x-direction, -i (u^ + Ug) 


dimensionless load parameter, f/e ' R 2 
Cartesian coordinates 
pressure-viscosity coefficient 
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Vr 


$ 


V* 

fluid viscosity at 1-atmosphere pressure 


Poisson's ratio 
v-1 


1 +■ 


30 


Subscripts: 

A solid A 

B solid B 

c central 

E elastic 

I isoviscous 

min minimum 

R rigid 

V viscous 

x,y in x and y directions, respectively 


DIMENSIONLESS GROUPING 

The dimensionless groupings used herein were also used by the authors in earlier 
publications (refs. 7 to 9) and are 

Dimensionless film thickness: 


H = 


R„ 


(1) 


where 


R x r Ax r Bx 


( 2 ) 


Dimensionless speed parameter: 


4 



where 



Dimensionless load parameter: 



Dimensionless material parameter: 


( 4 ) 

( 5 ) 


(6) 


G =■— — — « aE' (7) 

^iv, as 

Ellipticity parameter: 

k=- (8) 

b 

The conventional method of calculating the ellipticity parameter is to find a solu- 
tion to a transcendental equation that relates the ellipticity parameter and the elliptic 
integrals of the first and second kind to the geometry of the contacting solids. This is 
usually accomplished by some iterative numerical procedure. The following simple 
expression for the ellipticity parameter, which eliminates the necessity for that pro- 
cedure, is outlined in reference 10: 


k = 1 .03 



0.64 


( 9 ) 


where 



1 


( 10 ) 


*v 


1 


+ 


r Ay r By 


The approximate solution of the ellipticity parameter as obtained from equation (9) is 
within 3 percent of the exact solution for ellipticity parameters between 1 and 10 . 

The dimensionless group {H, U, W, G, k} has been a useful tool in understand- 
ing the results of references 7 to 9. Several authors (refs. 1 to 6) have noted that this 
dimensionless group can be reduced by one parameter -without any loss of generality - 
by using dimensional analysis. The film -thickness contours for the four fluid-film 
lubrication regimes are best represented graphically by the fewest parameters. 

Johnson (ref. 5) points out that the behavior distinguishing the four lubrication 
regimes can be characterized by three quantities, each having the dimensions of pres- 
sure, namely: 

(1) The reduced pressure parameter q^, a measure of the fluid pressure gener- 

ated by an isoviscous lubricant when elastic deformation is neglected 

(2) The inverse pressure- viscosity coefficient 1 /a, a measure of the change of 

viscosity with pressure 

(3) The maximum Hertzian pressure p^, the maximum pressure of a dry elastic 

contact 

Although Johnson's paper (ref. 5) does not consider elliptical contacts, it does state 
what the nondimensional parameters would be, namely: 


Dimensionless film parameter: 

-Or/ 


(11) 


Dimensionless viscosity parameter: 


_ • ^ _ GW 

>v = ac k T 

IT 


Dimensionless elasticity parameter: 


g E 


_ 1 /sV^ 3 /VV W 8 ^ 3 


ir \2> 


U 


( 12 ) 


(13) 


The ellipticity parameter k remains as discussed in equations (8) and (9). Therefore 

A, 

the reduced dimensionless group is {H, gy, g E , k} . 
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I so viscous-Rigid Regime 


The influence of geometry on the isothermal hydrodynamic film separating two 
rigid solids was investigated in reference 9 for fully flooded isoviscous conditions. 

The effect of geometry on the film thickness was determined by varying the radius 
ratio R /R from 1 (a circular configuration) to 36 (a configuration approaching a 

y x 

rectangular contact). The film thickness was varied over two orders of magnitude 
for steel solids separated by a paraffinic mineral oil. It was found that the minimum 
film thickness had the same speed, viscosity, and load dependence as the classical 
Kapitza (ref. 11) solution. However, incorporating the Reynolds cavitation boundary 
condition (i. e. , 3p/9n = 0 and p = 0 at the cavitation boundary, where n represents 
the normal coordinate to the cavitation boundary) resulted in an additional geometry ef- 
fect. Therefore, from reference 9 the dimensionless minimum (or central) film pa- 
rameter for the isoviscous-rigid lubrication regime can be written as 


{H . } ={H } 

L min-'iR L c j ir 


128/3$ 2 0.13 tan 


|+ 1.68 


where 


(14) 


,3 

R„ 


1.03 


1/0.64 


(15) 


$ = 



2_ 

3/3 



(16) 


In equation (14) the dimensionless film parameter is strictly a function of the geometry 
of the contact R/R x . 


Viscous-Rigid Regime 


From Blok (ref. 12) we can write the minimum film thickness for the viscous-rigid 
lubrication regime in a rectangular contact as 

\l/3 

^J!"t)"u"R 


h min = h c = 1.66(of 2 ^u 2 R x ) 1 


(17) 


By using the elliptical effect from reference 7 , equation (17) can be rewritten as 
h min = h c * 16 « d - e ‘°' 68k ) 


(18) 
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Note the absence of an applied-load term in equation (18). When expressed in terms of 
the dimensionless parameters of equations (11) and (12), this can be written as 


fW,,, = l«c> VR = 166 SV 3 * 1 - e '°' 68k > 


(19) 


Note the absence of the dimensionless elasticity parameter in equation (19). 


Isoviscous-Elastic Regime 


The influence of the ellipticity parameter k and the dimensionless speed U, load 
W, and material G parameters on the minimum and central film thicknesses was in- 
vestigated theoretically for the isoviscous-elastic regime in reference 8 . The elliptic- 
ity parameter was varied from 1 (a circular configuration) to 12 (a configuration ap- 
proaching a rectangular contact). The dimensionless speed and load parameters were 
each varied by one order of magnitude. Seventeen cases were considered in obtaining 
the dimensionless minimum-film-thickness equation 

H min = 7 - 43 U°- 65 W'°- 21 (1 - 0.85 e“°- 31k ) (20) 

From equations (11) and (13) the general form of the dimensionless minimum-film- 
parameter equation for the isoviscous-elastic lubrication regime can be expressed as 

^min = A 4( 1 -°- 85e "°' 31k ) ( 21 ) 


where A and c are constants to be determined. From equations (11) and (13), equa- 
tion (21) becomes 


H ■ = AU 2 - 2 V 8/3C >- 2 (1 - 0.85 e -0,31k ) 

mm v • 


( 22 ) 


Comparing equation (20) with equation (22) gives c = 0. 67 . Substituting this into equa- 
tion (21) while solving for A gives 


A = 


H 

min 


gg ,67 (l - 0.85 e °- 31k ) 


(23) 


The arithmetic mean for A based upon the 17 cases considered in reference 8 is 
8.70, with a standard deviation of ±9.05. Therefore the dimensionless minimum film 
parameter for the isoviscous-elastic lubrication regime can be written as 
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(24) 


{« min> IE ' 8 ' 70 4' ' 67 ( 1 - °' 85 e '°' 31k > 


With a similar approach the dimensionless central film parameter for the isoviscous- 
elastic lubrication regime can be written as 



11 . 15 g° * 67 (1 - 0 . 72 e"° • 28k ) 


(25) 


Viscous-Elastic Eegime 

In reference 7 the influence of the ellipticity parameter k and the dimensionless 
speed U, load W, and material G parameters on the minimum and central film 
thicknesses was investigated theoretically for the viscous-elastic regime. The ellip- 
ticity parameter was varied from 1 to 8, the dimensionless speed parameter was var- 
ied over nearly two orders of magnitude, and the dimensionless load parameter was 
varied over one order of magnitude. Conditions corresponding to the use of solid mate- 
rials of bronze, steel, and silicon nitride and lubricants of paraffinic and naphthenic 
oils were considered in obtaining the exponent on the dimensionless materials parame- 
ter. Thirty -four cases were used in obtaining the following dimensionless minimum- 
film -thickness formula 


H . = 3.63 u°* 68 G 0,49 W 0,073 (1 _ e "°- 68k ) 
min v 7 


(26) 


The general form of the dimensionless film-parameter equation for the viscous- 
elastic lubrication regime can be written as 




(27) 


where B, d, and f are constants to be determined. From equations (11), (12), and 
(13) equation (27) can be rewritten as 


H - = BG d U 2_2d 2f W 2+3d+ (8f/3) (1 _ e -0.68k 
mm v 7 


(28) 


Comparing equation (26) with equation (28) gives d = 0.49 and f = 0.17. Substituting 
these values into equation (27) and solving for B give 


B = 


H 


0.49 0.17 
gy gg 


min 
(1 - e 


-0.68k) 


(29) 
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For the 34 cases used in reference 7 in obtaining equation (26), the arithmetic mean 
for B was 3.42, with a standard derivation of ±0.03. Therefore, the dimensionless 
minimum film parameter for the viscous-elastic lubrication regime can be written as 




3.42 gy 49 g° ,17 (l - e °- 68k ) 


(30) 


An interesting observation to make in comparing equations (19), (24), and (30) is that 
the sum of the exponents on gy and gj, is close to the value of 2/3 required for com- 
plete dimensional representation of these three lubrication regimes: viscous-rigid, 
isoviscous-elastic, and viscous-elastic. 

With an approach like that just described, the dimensionless central film param- 
eter for the viscous-elastic lubrication regime can be written as 


{ H c> 


VE 


3.61 g 


0.53 0.13 
V Se 


(1 


0.61 


e-°- 73k ) 


(31) 


PROCEDURE FOR MAPPING LUBRICATION REGIMES 


Having expressed the dimensionless minimum film parameters for the four lubri- 
cation regimes in equations (14), (19), (24), and (30), we used these equations to de- 
velop a map of the lubrication regimes in the form of dimensionless film-parameter 
contours. These maps are shown in figures 2 to 6 on a log-log grid of the dimension- 
less viscosity and elasticity parameters for ellipticity parameters of 1, 2, 3, 4, and 6, 
respectively. The procedure used to obtain these figures was as follows: 

(1) For a given ellipticity parameter k, {H min } was calculated from equa- 

IR 


tion (14). 

(2) For a value of H •„ ^ {H 
' ' mm 1 r 

step 1, the dimensionless viscosity parameter gy was calculated from equation (19) 


Tninj-jR and the ellipticity parameter k chosen in 


as 



TJ 

min 


3/2 


[_1. 66(1 - e °- 68k )j 


(32) 


This established the dimensionless film-parameter contour H m ^ n in the viscous-rigid 
regime. 

(3) For the values of k selected in step 1, H ^ selected in step 2, and gy ob- 
tained from equation (32), the dimensionless elasticity parameter was calculated from 
the following equation, which was obtained from equation (30): 
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-|1/0.17 


g E 


H 


min 


| 3 . 42 gy 49 (1 - e ~°- 68k ) 


(33) 


This established the boundary between the viscous- rigid and viscous-elastic regimes. 

A 

(4) For the values of k and H. chosen in steps 1 and 2. the dimensionless 
elasticity parameter was calculated from the following equation obtained by rearrange- 
ment of equation (24): 


S E 


H min 


“1 1/0. 67 


8.70(1 - 0.85 e"°- 31k ) 


(34) 


This established the dimensionless film-parameter contour H . in the isoviscous- 

r mm 

elastic lubrication regime. 

(5) For the values of k and H min selected in steps 1 and 2 and the value of g E 
obtained from equation (34), gy was calculated from the following equation: 




H 


“1 1/0. 49 


mm 


)_3. 42 g 9 ’ 17 (1 - e -0 - 68k )_ 


(35) 


This established the isoviscous-elastic and viscous-elastic boundaries for the particu- 

A 

lar values of k and H, - chosen in steps 1 and 2. 

mm c a 

(6) At this point, for particular values k and H the contours were drawn 
through the viscous-rigid, viscous-elastic, and isoviscous-elastic regimes. A new 
H m in va l ue was then selected, and the contour was constructed by returning to step 2 . 
This procedure was continued until an adequate number of contours had been generated. 
A similar procedure was followed for the range of k values considered. 


DISCUSSION OF CONTOUR PLOTS 

The maps of the lubrication regimes shown in figures 2 to 6 were generated by 
following the procedure outlined in the previous section. The dimensionless-film- 

A 

parameter H ^ contours were plotted on a log-log grid of the dimensionless vis- 
cosity parameter gy and the dimensionless elasticity parameter g £ for ellipticity 
parameters of 1, 2, 3, 4, and 6, respectively. The four lubrication regimes are 
clearly shown in these figures. The smallest film-parameter contour represents the 
value obtained from equation (14) and is the iso viscous- rigid boundary. The value of 
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this dimensionless film parameter for the iso viscous- rigid boundary increases as the 
ellipticity parameter k increases. 

By using figures 2 to 6 for given values of the ellipticity parameter k, the vis- 
cosity parameter gy, and the elasticity parameter g E , the lubrication regime can be 
ascertained and the approximate value of the dimensionless film parameter H ^ de- 
termined. When the lubrication regime is known, a more accurate value of H m ^ n can 
be obtained by using the appropriate dimensionless-film-parameter equation. 

Figures 7 to 11 were obtained by using the dimensionless-film-parameter equa- 
tions for the four lubrication regimes and a procedure like that used in obtaining fig- 
ures 2 to 6. The difference between the two groups of figures is that the ordinate in 
figures 2 to 6 is the viscosity parameter g v and the contours are values of the dimen- 

A ' A 

sionless film parameter H m ^ n but the ordinate in figures 7 to 11 is H min and the 
contours are values of gg . In figures 7 to 11 the boundaries of the isoviscous-rigid 
and isoviscous-elastic regimes are shown, as are the viscous-rigid and viscous- 
elastic regimes. 

Figure 12 gives a three-dimensional view of the surfaces of three constant dimen- 
sionless film parameters = 500, 2000, and 6000). The coordinates in this fig- 

ure are the dimensionless elasticity parameter g E , the dimensionless viscosity pa- 
rameter gy, and the ellipticity parameter k. The four fluid-film lubrication regimes 
are clearly shown. In this figure, we attempt to show how the parameters gy, g E , 
and k influence the dimensionless film parameter H min as well as to define the 
fluid-film lubrication regimes. 


CONCLUSIONS 

The dimensionless minimum-film-parameter H min equations for the four lubri- 
cation regimes found in elliptical contacts are 


Isoviscous-rigid regime: 


where 


{H . } = 128/3$ 

mm IR 


.2 


0 . 13 tan' 





1/0.64 


12 



Viscous-rigid regime: 


-ter 


{“mini™ ' 1 -«« «v V - e-°' 68k ) 


VR 


Isoviscous-elastic regime; 




Viscous-elastic regime: 


{»min> 


0 . 490 . 17 , 


VE 


= 3.45 gy gg (1 - 


-0.68k 


where and Ry are the effective radii in the x and y directions, k is the ellip- 
ticity parameter, gy is the dimensionless viscosity parameter, and g^, is the dimen- 
sionless elasticity parameter. The relative importance of the elastic and viscous ef- 
fects is the factor that distinguishes these regimes. 

In addition these equations have been used to map the lubrication regimes by plot- 
ting film-thickness contours on a log-log grid of the viscosity and elasticity parameters 
for five values of the ellipticity parameter. These results present, for the first time, 
a complete theoretical film-parameter solution for elliptical contacts in the four lubri- 
cation regimes. The results are expected to be most useful in the solution of many 
practical lubrication problems. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, June 29, 1978, 

505-04. 
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Figure 2. - Map of lubrication regimes with dimensionless-film-parameter contours on log-log grid of 
dimensionless viscosity and elasticity parameters for elliplicity parameter of l(k » 1). 
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Figure 5. - Map of lubrication regimes with di men si on less-fi I m-parameter contours on log-log grid of 
dimensionless viscosity and elasticity parameters for ellipticity parameter of 4 (k - 4). 



Dimensionless viscosity parameter, 







Dimensionless minimum film parameter, H m j n Dimensionless minimum film parameter, 



Dimensionless elasticity parameter, g^ 

Figure 7. - Map of lubrication regimes with dimensionless-viscosity-parameter contours on log-log grid of 
dimensionless film and elasticity parameters for ellipticity parameter of 1 (k = 1). 



Figure 8. - Map of lubrication regimes with dimensionless-viscosity-parameter contours on log-log grid of 
dimensionless film and elasticity parameters for ellipticity parameter of 2 (k = 2 ). 



Dimensionless minimum film parameter, H min Dimensionless minimum film parameter, H 



10 2 10 3 10 * 10 5 10 6 10 7 
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Figure 9. - Map of lubrication regimes with dimensionless-viscosity-parameter contours on log-log grid of 
dimensionless film and elasticity parameters for ellipticity parameter of 3 (k - 3). 



Figure 10. - Map of lubrication regimes with dimensionless-viscosity-parameter contours on log-log grid of 
dimensionless film and elasticity parameters for elliplicity parameter of 4 (k - 4). 
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Figure 11. - Map of lubrication regimes with dimensionless-viscosity-parameter contours on log-log grid of 
dimensionless film and elasticity parameters for ellipticity parameter of 6 (k = 6). 
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